Wigner crystal induced by dipole-dipole interaction in one-dimensional optical lattices 
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We demonstrate that the static structure factor, momentum distribution and density distribu- 
tion provide clear signatures of the emergence of Wigner crystal for the fermionic dipolar gas with 
strongly repulsive dipole-dipole interactions trapped in one-dimensional optical lattices. Our numer- 
ical evidences are based on the exact diagonalization of the microscopic effective lattice Hamiltonian 
of few particles interacting with long-range interactions. As a comparison, we also study the system 
with only nearest-neighbor interactions, which displays quite different behaviors from the dipolar 
system in the regime of strong repulsion. 
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I. INTRODUCTION 

Ultracold dipolar atomic and molecular gases with 
long-range interactions have become a very active re- 
search field of cold atom physics in past years. Dipolar 
interactions with long-range anisotropic character have 
been observed in Chromium atoms [l|, 0]. More re- 
cently, the development of stimulated Raman adiabatic 
passage technique has succeeded in creating a nearly de- 
generate gas of heteronuclear polar molecule gases which 
have much greater dipole momentum and thus are 
promising candidates with very strong dipolar interac- 
tions. The experimental progress has stimulated theoret- 
ical studies of dipolar gases with long-range anisotropic 
interactions. Loading ultracold polar molecules into opti- 
cal lattices also provides a fascinating platform for study- 
ing quantum many-body systems with long range inter- 
actions. 

In comparison with the short-range and isotropic in- 
teractions, dipole-dipole interactions (DDIs) may induce 
many new effects and phases in quantum gases due to 
their long-range and anisotropic character. There have 
been many theoretical works on, for example, the po- 
larized dipolar Bose-Einstein condensations (BECs) [g- 
8) , spinor-dipolar BECs foj-TTlT] , supersolid (l3jl6j , s- wave 
scattering resonances [17fl, quantum Hall effects [18| and 
one-dimensional (ID) dipolar systems [l9l - l23j . In the 
presence of long range interaction, an interesting issue 
is the emergence of the Wigner crystal phase. A Wigner 
crystal is a crystalline phase of electrons first predicted by 
Wigner [24j]. As the interacting potential energy domi- 
nates the kinetic energy at low densities, the particles 
tend to form a regular crystal due to the strong long- 
range interaction. Wigner crystal in one dimension with 
long range interactions has been theoretically studie d by 
using different methods, including bosonization |25l - t28| , 
quantum Monte Carlo (QMC) [Io|, [H US HI and exact 
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diagonaliztion 22] . In the seminal work of Schulz |25| , it 
has been clarified that the Coulomb long-range repulsion 
enhances the 4kp charge density correlations and drives 
the system to the Wigner crystal phase. 

In this paper, we study the fermionic dipolar gas with 
strong dipole-dipole interactions trapped in ID deep op- 
tical lattices, which can be effectively described by a 
Fermi lattice model with long-range interactions falling 
off as 1/x 3 . For ID fermionic models, the case of interac- 
tions falling off as l/x a has been studied in the scheme 
of bosonization (2(| l27j . Treating long-range forward 
scattering as a perturbation, they find that for a > 1 
the long-range forward scattering is an irrelevant per- 
turbation. As the above conclusion is obtained based 
on the perturbation analysis of the low-energy effective 
Luttinger liquid theory with linear dispersion, it does not 
exclude the existence of Wigner crystal phase in the pres- 
ence of strong dipole-dipole interactions. In this work, we 
shall study the ground state properties of the Fermi lat- 
tice model with long-range interactions by means of the 
exact diagonalization method. To see clearly the effect 
of the long-range interaction, we carefully compare the 
long-range interacting dipolar systems against the Fermi 
systems with short-range interactions. Our results dis- 
play significantly different behaviors for the Fermi sys- 
tems with long-range and short-range interactions in the 
strongly interacting regime. By observing the differences 
of static structure factor, momentum distribution and 
density distribution, we give clear evidences for the emer- 
gence of Wigner crystal for dipolar fermions with strong 
dipole-dipole interactions. 

The paper is organized as follows. In section II, we de- 
rive the effective lattice model for the fermionic dipolar 
gas in a ID optical lattice. In section III, we present our 
results for both the dipolar Fermi system with long range 
interaction and system with only nearest-neighbor (NN) 
interaction by using the exact diagonalization method 
and compare their different behaviors with increasing the 
interaction strength. A summary is given in the last sec- 
tion. 
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II. MODEL HAMILTONIAN 



Let us consider a system of electric or magnetic dipoles 
confined in a quasi-lD optical lattice along the x direc- 
tion. The interaction between dipoles aligned in the x— z 
plane by a homogeneous hied can be simplified to 



V dd (r) = C, 



1 - 3 cos 2 6L 



dd- 



(1) 



where cos 9 rc i = r ■ d/(rd) and Cdd measures the strength 
of the DDI. The strength of the DDI is given by 
Cdd = rf 2 /(47T£ ) for two electric dipoles, and by Cdd = 
/io<i 2 /(47r) for the magnetic ones, with Eq and fiQ be- 
ing the electric constant and magnetic constant, re- 
spectively. In the single-mode approximation with only 
the transversal ground state considered, we can inte- 
grate over transversal directions and get the effective 
ID DDI given by V dd = U dd V dd (\x\/l ± ), with U dd = 
-C dd [l + 3cos(20)]/(8Zi), V dd {u) = -2u + v^F(l + 
u 2 ) exp(w 2 /2)erfc(w/v / 2) and l± = y/h/(rruj±) [13, Hi- 

The Hamiltonian for dipolar fermions in optical lattices 
is given by 



II = I dx^(x) 
1 



_&i_&_ + v(x] 
2m dx 2 ' 



4,{x) (2) 
dxdx'^(x)ip^(x')V dd (x - x')ip(x')ij}(x), 



where m is the mass of the dipolar fermions, ip{x) 
is a fermionic field operator for the dipolar fermion, 
and Vq(x) is the optical lattice field given by Vq(x) = 
Vosin 2 (fca;) with the wavevectors k = 2ir/\ and the 
wavelength of the laser light A. Considering the deep 
optical lattice with particles trapped in the lowest vi- 
brational state lo{x) = exp(— a; 2 /(2Z 2 ))/(7r 1 / 4 v / U) with 
l x = y/H/ (mujx) and ftw x = 2^/ErVq, we use the Wannier 
basis to expand the field operator ip(x) = Ciiv(x — Xi). 
Similar to the case of Hubbard model [31] , we can get the 
effective lattice Hamiltonian for the polar fermionic sys- 
tem 

H = - J^>]c i+1 + H.c.) + 1 ]T V dd {\i - j\)n inj , (3) 

where cj(ci) is the creation (annihilation) operator of the 
fermion, J = s/ErVq cxp(— n 2 \/Vq/E r /4) with E R — 
H 2 k 2 /2m the recoil energy of the system, — j\) = 

Udd J J &x&x'V dd {\x-x'\/lx)bJ*{x-Xi)w*{x' ' — Xj)u)(x — 
Xi)u>(x' — Xj) is the DDI between the dipoles at the po- 
sitions of Xi and Xj. 

After the numerical calculation by setting Er, h, and 
m as units, lj±/uj x = a and Vq/Er = j3, the V dd (\i—j\) is 
found to decay as l/\i — j\ 3 according to FigfTJ In Fig. [U 
we show the Log-Log plot of Vdd/v d d{0) versus distances 
between dipoles in ID optical lattices for different a and 
/3 with v d d(0) = — C d d[l + 3 cos(20)], where 9 is the angle 
between the dipole direction and x axis. When a is fixed 
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FIG. 1: (Color online) Log-Log plot of Vdd/vdd(6) versus dis- 
tances between two dipoles in optical lattices for different a 
and f3 with v dd (0) = ~C d d[l + 3cos(20)]. (a) We set a = 20 
for different /3, and the slop of the solid line is —3; (b) /3 = 20 
is fixed for a = 15, 20, 25, also the gradient of the fitting line 
is —3. 



for different (3 (see Figflja)), the gradient of the fitting 
line is —3. In the Figfljb), the slops of the fitting line is 
the same as the one in the FigOJa) where (3 is fixed and 
a is changing. Also we can see when a is large enough, 
Vdd changes little as a increases, because all the dipoles 
have been confined along the axial direction and can be 
seen as a one dimensional system. So we can simplify the 
Hamiltonian ([3]) as 



H= -J^(cjc l+1 +H.c.) 



2 ^ i 



J\ 



(4) 



with V the strength of the DDI. For comparison, we also 
consider the system with only short-range interactions 
described by the following hamiltonian 



Xh% TLj j 



(5) 



where only the nearest-neighbor interaction is considered 
and means summation over nearest neighbors. In 

the present work, we only consider the case with repulsive 
interaction V > 0. For simplicity, we take 6 = n/2. 



III. RESULTS 

We shall study the ground state properties for both 
systems described by Hamiltonians (0} and ([5]) by the ex- 
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FIG. 2: (Color online) Static structure factors S(k) for sys- 
tems with L = 30, N — 5 and different V. (a) The dipolar 
Fermi system, and (b) the Fermi system with only the nearest- 
neighbor interaction. 



FIG. 3: (Color online) Static structure factors S(k) vs k/kF 
for systems with L — 30, V = 5000 and different N. (a) is 
for the dipolar system, and (b) is for system with short range 
interaction. 



act diagonalization method. To give a concrete example, 
we focus on the system with five particles (N = 5) in a 
lattice with size L — 30 under the periodic boundary con- 
dition. In order to gain some intuitive understanding of 
effect of long-range interactions, we first consider the lim- 
iting cases before presenting our calculated results. In the 
limit of V/ J — > 0, the hopping term dominates and pre- 
vents the formation of crystal phase, and thus the differ- 
ence of the long-range interaction and nearest-neighbor 
interaction is not obvious. However, in the strongly in- 
teracting limit of V/J — > oo, the hopping term can be 
ignored, and the effect of long-range interaction becomes 
significant. The long-range repulsive interaction tends to 
repel particles to be equally spaced and form a solid like 
state. When apart from the strongly interacting limit, 
the hopping processes and quantum fluctuations prevent 
the formation of a perfect Wigner crystal. 

To characterize the phase of Wigner crystal, we calcu- 
late the static structure factor, defined as 

S ( k ) = \Y.j k{i - j) [{mn j ) ~ (rnXn,-)], (6) 

where k = 2rmr/L with m = 0,1,..., L. The phase of 
Wigner crystal can be characterized by the onset of the 
Akp peak in the static structure factor with kp = mr 
and n = N/L. When V is small, we do not find obvi- 
ous differences of S(k) for systems with either long-range 
or short-range interaction. However, significant differ- 
ences are emergent as the interaction strength increases 



to the strongly interacting limit. In Figure [2] we dis- 
play the static structure factor S(k) for different values 
of V = 0, 10, 50, 100 and 500. For both systems, there is a 
peak at k = 2kp for V — 10 despite that the peak for the 
system with long-range interaction is more obvious. As 
the interaction increases further, more peaks emerge at 
reciprocal lattice vectors kL/2irN = integer for the sys- 
tem with the long-range interaction (Figj2ji) . The height 
of peak increases with increasing V, evidencing the occur- 
rence of the Wigner crystal state. While for the system 
with short range interaction described by Hamiltonian 
([5]), from FigJSJb) we can see that only peaks at ±2kp 
occur and there are no other peaks appearing. Also, the 
height of peaks at ±2kp increases more slowly with the 
increase of V in contrast to the system with dipole-dipole 
interaction. 

We also demonstrate data of static structure factor 
S(k) versus k/kF in Fig[3] for systems with L = 30, 
V = 5000 and N = 2,3,4,5. Systems with different N 
display similar behaviors in the large V limit. As shown 
in Figj3ja), for systems described by Eq.Q with different 
filling factors, peaks of S(k) emerge at reciprocal lattice 
vectors k/kF = 2m with m the integer. The appearance 
of peaks for various filling cases indicates that the Wigner 
crystal emerges when the filling is either commensurate 
or incommensurate with the optical lattice. For the sys- 
tems described by Eq.©, FigJ^b) shows that there are 
no other peaks except at k/kp = 2. 

Next we calculate the momentum distribution, which 
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FIG. 4: (Color online) Momentum distributions for systems 
with L = 30, N = 5 and different V. (a) The dipolar Fermi 
system, (b) the Fermi system with nearest-neighbor interac- 
tion. 
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FIG. 5: (Color online) Density distributions for systems with 
L = 30, N — 5 and different V under the open boundary 
condition, (a) The dipolar Fermi system, and (b) the Fermi 
system with nearest-neighbor interaction. 



is defined by the Fourier transform with respect to i — j 
of the one-particle density matrix with the form 



(7) 



In FigHl we show momentum distributions for both sys- 



tems with different V. Fig|4}a) is corresponding to the 
dipolar model, and FigHJb) is for the model with NN in- 
teractions. We also find the similarity of the two different 
models in the regime of small V. They show typical mo- 
mentum step distribution of fermions. But with increas- 
ing V , the momentum densities of the dipolar model be- 
come fatter and lower (FigHfa)). At last for very strong 
V, the momentum distribution displays the Gaussian dis- 
tribution which mainly comes from the Fourier transfor- 
mation of the diagonal part of the density matrix, reveal- 
ing that there is no overlap between the localized wave 
packets with the non-diagonal density matrix tending to 
vanish. In contrast to the case with DDI, the momen- 
tum distributions in Fig@Jb) change not obviously even 
for very strong V. 

To gain an intuitive insight on the crystalline phase, we 
analyze the density distribution m = (c\ci) for systems 
described by Hamiltonian ^ and ([5]) with open bound- 
ary conditions. In Fig. [5] we show density distributions 
of five particles for different V. As shown in FigfSta) 
for the DDI systems, with increasing the long-range in- 
teraction, particles move apart each other with more 
sharp peaks emerging in the density profile. When V 
is large enough, the density resembles five well-separated 
localized wave packs, which are totally apart from each 
other with the equilibrium positions of wave packets min- 
imizing the interaction energy, characterizing the system 
evolved into the Wigner crystal regime [22j . As a com- 
parison, the density distributions shown in FigJSJb) only 
change marginally with increasing the nearest-neighbor 
repulsion. No crystalline signal is detected even in the 
limit of strong repulsion. 

Before ending this paper, we would like to discuss the 
validity of the approximation of the lowest transverse 
state in deriving the effective ID Hamiltonian @. Ex- 
cept of the requirement of lo±/lu x ^> 1, one also needs 
that the interacting energy Eddi induced by the dipo- 
lar interaction is much smaller than the energy of trans- 
verse confinement fkj±. Different from the case with 
contact interaction, the interacting energy of the dipo- 
lar system is proportional to the interaction strength V, 
and thus the approximation is expected to break down 
if Eddi > fujj±. For the present case, the longitudinal 
kinetic energy is greatly suppressed due to the existence 
of the longitudinal optical lattice, i.e., J -C tko x . Con- 
sequently, even when V/J 3> 1, the interacting energy 
is still much smaller than Hoj± and thus the effective ID 
description still holds true even in the regime of forma- 
tion of Wigner crystal. Taking the case corresponding 
to Fig. 2 as an example, if we take a = 20 and j3 = 20, 
J/huj x ~ 10~ 6 , and therefore even for V/J = 500 the 
Eddi is still much smaller than hu±. For systems with 
even stronger interaction or with large atom numbers, 
one can also tune a > 1 to fulfill the requirement of 
Eddi "C huj±. In the future work, it would be also inter- 
esting to study the crossover from ID to high-dimensional 
system when the transverse confinement decreases, for 
which the approximation of the lowest transverse state 
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does not work well and new phenomena may appear [32j ] . 

IV. SUMMARY 

In summary, we have studied the dipolar fermionic sys- 
tem with dipole-dipole interactions trapped in ID optical 
lattices by means of the exact diagonalization method. 
We have shown that the static structure factor, momen- 
tum distribution and density distribution provide clear 
evidences for the existence of Wigner crystal, as the in- 
teraction energy overcomes the energy scale of hopping 
energy. We also compare our results of the dipolar sys- 
tem to the model with only nearest neighbor interactions, 



which exhibits no signatures of Wigner crystal even in the 
strongly interacting limit. Our study unveils the impor- 
tant role of the long range interaction in the formation 
of the Wigner crystal. 
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